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Abstract —In the preliminary design process of switched 
reluctance machine (SRM), adjusting the number of wingding 
turns or cross-sectional area of conductor to optimize windings 
design has dramatic effects on the performances like torque 
density, efficiency and thermal dissipation and so on. However, 
the difficulty exists on how to guarantee these performances 
while optimizing winding design. This paper proposes an 
analytical optimization design method based on an regulation 
model of ampere density and coil space factor, which can directly 
determine the optimal number of winding turns, the cross- 
sectional area of conductor, the ampere density and the coil space 
factor. And a MATLAB pre-design program has been developed 
to provide two design schemes respectively with and without the 
optimization model. The comparison analysis has been further 
carried out with finite element analysis (FEA). The simulation 
results verify that the performances can be highly improved with 
the proposed optimization model. 
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ampere density, coil space factor Introduction 


I. INTRODUCTION 


Switched reluctance machine (SRM) has the concentrated 
windings on the stator but no permanent magnets and 
windings on the rotor. It is an attractive alternative for electric 
vehicle (EV) and hybrid electric vehicle (HEV) applications 
owing to its simple and low cost construction, wide constant 
power range, low inertia and low inrush current et al[1]. 

In many applications, it demands for high torque density, 
high efficiency, high output power, excellent thermal 
dissipation and low cost, which can be improved through 
comprehensive optimization design of SRM and its optimal 
control. Because the excitation and unique inductance profile 
is closely determined by the rotor position, the mathematical 
model of SRM is highly nonlinear and the flux linkage is easy 
to saturated at the poles, which makes it harder to design and 
adjust the geometric dimensions to meet the design 
specifications. Some works focus on the comprehensive 
design procedure and program for common SRM design, 
which allows an efficient and fast design by understanding the 
influence of relevant design parameters. Berker Bilgin et al. 
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have proposed a MATLAB based design tool to calculate the 
geometry with the conventional magnetic circuit method and 
automatically draws and simulates the SRM model in FEA 
software without considering the optimization design[2]. And 
some works focus on the geometry and configuration 
optimization to improve the performances like efficiency, 
torque density and power output et al. Kazuhiro Ohyama et al. 
have studied the approach to improve the efficiency through 
designing the cross sections and axial shapes of rotor and 
stator cores using the magnetic field analysis with 2D and 3D 
finite element methods[3]. J. Hur et al. have implemented the 
flux barriers to improve the torque density[4]. Anas Labak et 
al. have proposed a design and analysis method of a novel 
axial-flux SRM, in which the detailed procedures of deriving 
the output power equation as a function of the motor 
dimensions and parameters and an exclusive pole-shape 
design are presented. It focuses on the structural design to 
improve the performance[5]. Dong-Hee Lee et al. have 
proposed the first and second optimization process to reduce 
the torque ripple and improve the output torque through 
optimizing the air gap which is considered to be a function of 
the rotor position[6]. 

However, the empirical values are generally needed to 
guide and amend the design process which always needs 
abundant design experiences. One of the crucial design 
procedure is to design and optimize the winding, i.e., the 
number of winding turns, slot fill factor, ampere density and 
wire diameter. And winding design and optimization also has 
great effects on the SRM performances[7-9]. Therefore, this 
paper mainly focuses on the winding design and optimization 
issue, and proposes a optimization regulation model based on 
the ampere density and coil space factor optimization, which 
is actually a winding design and optimization approach, 
aiming at improving the torque density and the output power 
as well as the efficiency through choosing optimal number of 
winding turns, wire diameter, ampere density and coil space 
factor. An comparison study with magnetic circuit method and 
FEM simulation also has been carried out. It verifies that the 
proposed optimization model is not only advantageous for 
improvement of torque density, but also for power output and 
efficiency. 


Il. OPTIMIZATION REGULATION MODEL 


Ampere density and coil space factor are two crucial 
factors that must be comprehensively investigated and 
designed in certain appropriate intervals. It means that , on the 
other hand, the ampere density and coil space factor are non- 
unique. There are various combinations of these two factors to 
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meet the same design specifications and requirements. 
Therefore, guessing that it exists an optimal point, the problem 
is simply how to determine it. However, actually, whether 
there is an optimal point or not, it depends on the optimizable 
range, which is further determined by the preliminary 
structural design of SRM. The ampere density and coil space 
factor can be expressed as 
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where J is the ampere density, ks is the coil space factor, Npa is 
the windings turns in series of each phase, Da is the rotor 
outside diameter, A is the current per unit length, q is the 
number of phase, / is the effective winding current, Sw is the 
stator yoke window area, Seu is the net conductor cross- 
sectional area of each slot and S, is the cross area of conductor. 
Corresponding definitions are illustrated in Figure 1. And the 
relationship between the ampere density and coil space factor 
is shown as Figure 2. 
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Figure 1 Definitions of ampere density and coil space factor 
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Figure 2 Relationship between ampere density and coil space factor 
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From (1) and (2), it is obvious that the ampere density and 
coil space factor can be regulated by changing the number of 
winding turns and cross-sectional area of conductor assuming 
that Da, A, q and Sw change little and can approximately be 
considered to be constant. It hence yields that 
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Only when the intervals of J and ks having overlapping 
section on the curve, there is an optimal point. And the 
overlapping intervals is the so called optimizable range. More 
detail discussion of the optimizable range can be found in 
reference[10]. Once the optimization possibility has been 
discussed, another problem is how to find the optimal point. 
So, an optimized regulation model has been proposed, which 
adapts the electromagnetic torque Tem as the optimization 
objective function. According to the co-energy theory[11], the 
electromagnetic torque can be calculated with co-energy W”. 
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where N, is the number of rotor poles. Figures 3 shows the 
magnetization curves at aligned position @, and unaligned 
position @, , in which the shaded area represents the co-energy, 
where J, is the ideal square wave peak current, i, is the 
critical saturation current for y, (0, i), L, is the inductance 


at@ position, Lis the inductance at @, position when i< i. 


Figure 3 Magnetization curves at aligned (0, ) and unaligned (8, ) positions. 
Shown as Figure 3, the co-energy can be expressed as 
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On the other hand, the co-energy can also be expressed as a 
function of the number of winding turns in series of each 
phase [10, 12], 
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where P;~Ps are constants which are related to the flux and 
inductance at special rotor positions[10]. Thus, ultimately, 
the electromagnetic torque Tem can be expressed in the form 
of Npn as follows, 
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The practical task of the model is to determine the optimal 
ampere density J and coil space factor ks by simply regulating 
Sa and Nph in the optimization range, meanwhile guaranteeing 
providing optimal electromagnetic torque. Therefore, the 
optimized regulation model can be described as 
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There are also some constraint conditions as follows, 
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Thus, the optimization regulation model can ultimately 
be expressed as a standard nonlinear optimization model 
shown as 
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This is a standard nonlinear optimization model, which 
can be solved with various numerical computing methods like 
the Kuhn-Tucker (K-T) equation method and the Sequential 
Quadratic Programming (SQP) methods. There are also 
corresponding commands in MATLAB. 


HI. OPTIMIZATION DESIGN PROGRAM WITH MATLAB 


In order to verify the feasibility of the optimization 
regulation model, a pre-design program has been completed 
with MATLAB, which utilizes the conventional design 
approach and procedure to determine the initial geometric 
dimensions[13, 14]. It provides two preliminary designs, one 
of which is with the optimization regulation model, and the 
other is without the model. Figure 4 shows the program flow 
chart with the optimization regulation model. And the one 
with no optimization is the same with the procedure shown as 
in the dashed box. 

Corresponding preliminary design specifications and 
parameters without the optimization regulation model are 
listed in TABLE I. Compared with the design without the 


optimization 
parameters with the model are the same excepting the winding 
and electromagnetic parameters as shown in TABLE II. 


regulation model, the designed 


Initilize control parameters: 
turn-on angle, turn-off angle, 
conduction angle 


== as 


Caculate winding turns 


a Ss 


Modeling magnetization curves , 
caculate inductances , 
flux linkages 


Input 
design 
indicators 


Caculate rated voltage, 
electromagnetic power, 
torque 


Optimize 

turn-on & 

conduction] 
angles 


Caculate actual peak current, 
effective current, 
ideal square-wave current peak 


Initilize 
electromagnetic load 


Caculation winding parameters 


Caculate stator and 
rotor dimensions 


Caculate current density, 
coil space factor 


Caculate co-energy and designed 
electromagnetic torque 


Optimize rotor 
outer diameter 


Caculate air gap flux 
denstiy 


Designed electromangetic torque 
OK? 


Does air gap flux 
density fit demand 


Output optimization solution 


Figure 4 MATLAB program flow chart 


TABLE I Design specifications and parameters 


structural 


Design Specifications 


Number of phases 4 ee 8/6 
Rated voltage 380VAC Rated speed 1500rpm 
Rated output power 11KW Rated efficiency 88% 
Design Parameters 

Turn-on angle -2.35deg Turn-off angle 18.5deg 
Rotor inner diameter 59.5mm iee a 210mm 
Rotor yoke height 16.60mm aa 178.5mm 
Stator pole width 21.83mm_ | Stator pole arc 2ldeg 
Rotor pole width 23.72mm_ | Rotor pole arc 23deg 
Air gap length 0.4mm oo 119mm 
Stator slot depth 30.19mm | Stator yoke height | 14.19mm 
Stator pole pitch 62.31mm 


TABLE II Theoretical performances comparison 


Effects of the optimization regulation model 


Nominal Before Optimized | After Optimized 
Winding turns in series of 124 116 

each phase 

Cross area of 

conductor(mm2) aes 4:2087 
Wire diameter(mm) 2.24 25 
Ampere density((A/mm2) 5.53 4.79 

Coil space factor 0.5277 0.6051 
Rated electromagnetic 75.02 99.19 
torque(N.m) 
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The optimizable range is shown in Figure 5. The coil 
space factor ks and the ampere density J are expected to be 
designed in an certain appropriate interval, which is 
k, €[0.5,0.7] and J €[3.5,5.5] . After optimized, the coil 


space factor is k,;=0.6051, and the ampere density is 
J=4.79A/mm*. It means that both ks and J are in the 
optimizable range. And Figure 6 shows the theoretical 
relationship between the electromagnetic torque and the 
number of winding turns. It states that the number of winding 
turns in series of each phase N,, should be in the interval of 
N „ €[116,156] when the cross-sectional area of conductor is 


2=4.9087mm”, and the corresponding electromagnetic torque 
theoretically varies from 75.02N.m to 100.5N.m as shown in 
Figure 6. It means that the theoretical value of optimal 
electromagnetic torque should between 75.02N.m and 
100.5N.m, and TABLE II shows that the theoretical optimal 


electromagnetic torque is 99.19N.m, which is in the 
optimizable range. 
8+- . + 
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Figure 5 Ampere density vs. coil space factor (Optimizable range) 
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Figure 6 Electromagnetic torque vs. the number of winding turns 


It states that the optimization regulation model functions 
well. It can effectively search an theoretical optimal point in 
the optimizable range. Before optimized, the designed rated 
electromagnetic torque is theoretically about 75.02N.m, while 
after optimized, it increases to 99.19N.m, which can 
theoretically improved by 32.22%. 


IV. COMPARISON ANALYSIS 


And a performances comparison analysis has been carried 
out with ANSYS Maxwell. The purpose is merely 
comparatively study the effects of the optimization regulation 
model on the performances of SRM. The winding 
configuration here is the conventional type of long magnetic 
circuit loop shown as Figure 7, of which the corresponding 
excitation sequence A-A&B-B-B&C-C-C&D-D-D&A. 

Figure 8 shows the finite element mesh result, in which the 
maximum elements length of the coils is 9.02mm, and 9.1mm 
of stator and rotor. The maximum surface deviation is 
0.105mm, and the deviation angle is 15deg. 


(a) Phase A (b) Phase A&B (c) PhaseB (d) Phase B&C 
i E Q © © © 
Ce) PhaseC (f) PhaseC&D (g) PhaseD (h) Phase D&A 


Figure 7 Winding configuration and excitation phase sequence 


Y 


0 100 
Figure 8 Finite element mesh 


200 (mm) 


Figure 9(a)~(c) and Figure 10(a)~(c) respectively show the 
flux linkage and density distribution at certain special rotor 
positions including the partial overlapped, aligned and 
unaligned positions with respect to the pole of phase-A. 
Specifically, when the stator pole and the rotor overlapping, 
the flux density at corresponding overlapped part is relatively 
larger than that at aligned or unaligned positions. And TABLE 
II lists the average flux linkage and flux density comparison. 


TABLE III Flux linkage and flux density comparison 


Nominal Before Optimized After Optimized 
Flux linkage(Wb) 0.682204 0.640847 
Stator-pole flux density(T) 1.64944 1.65630 
Stator-yoke flux density(T) 1.14319 1.14795 
Rotor-pole flux density(T) 1.63482 1.64162 
Rotor-yoke flux density(T) 0.977221 0.981288 
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Figure 9 Flux linkage distribution(1500rpm) 
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Figure 10 Flux density distribution(1500rpm) 


TABLE IV lists the theoretical performances comparison. 
It shows that the performances has indeed improved after 
optimized. 
TABLE IV Simulation performances comparison 


Performances comparison 
Nominal Before optimized | After optimized 

Rated speed(rpm) 1761 1884 
Rated output torque(N.m) 69.35 69.97 
Output torque at 
1500rpm(N.m) 91.87 103.06 
Estimated start torque(N.m) 540.85 669.45 
Rated electromagnetic 
power(KW) 13.86 14.84 
Rated output power(KW) 12.79 13.81 
Maximum output 
power(KW) 16.91 18.41 
Winding copper loss(W) 688.10 597.01 
Total loss(W) 1074.55 1026.65 
Efficiency(%) 92.25 93.08 


Figure 11 shows the simulation output torque curve at 
different speed, which states that the output torque has been 
highly improved. When applied the same rated load toque of 
70N.m, it is obvious that after optimized, the steady speed 
increases from 1761rpm to 1884rpm, which means that the 
operation range has been improved with the optimization 
regulation model. Correspondingly, the rated output 
mechanical torque increased from 69.35N.m to 69.97N.m. 
Besides, before optimized, the output torque is 91.87N.m at 


1500rpm, while after optimized, it increases to 103.06N.m, 
which is improved by 12.18%. 

The simulation output power also has been improved 
shown as Figure 12. The output power at 1500rpm is about 
14.34KW before optimized and 16.09KW after optimized, 
which has been increased by 12.2%. 

Shown as Figure 13, it states that under the speed of 
1500rpm, the efficiency is not significantly improved, while 
above the rated speed it obviously increases compared with 
the circumstance before optimized. 
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Figure 11 Output torque comparison with various speed 
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Figure 12 Output power comparison with various speed 
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Figure 13 Efficiency comparison with various speed 


V. CONCLUSION 


This paper proposes a regulation model of ampere density 
and coil space factor to optimize the winding design of SRM. 
The purpose is to improve the performances through choosing 
the optimal number of winding turns, the cross-sectional area 
of conductor, which in turn leads to the optimal ampere 
density and coil space factor. It ultimately verifies that, after 
optimized, the ampere density is effectively decreased, which 
is beneficial for reducing the windings loss and increasing the 
efficiency. And the coil space factor is appropriately increased, 
which leads to the improvement of thermal performance of 
SRM on the other hand. Besides, not only the torque density, 
but also the power density and efficiency performances can be 
further improved. The simulation indicates that the average 
output power can be improved by 12.2% with the proposed 
optimization model and the rated efficiency increases from 
92.25% to 93.08%. 

The future work is torque ripple optimization and noise 
suppression through optimal design of SRM, and thermal 
performance is still remaining to be analyzed. 
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